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Abstract: We model the effect of photon and ultra-light pseudoscalar mixing on the propagation 
of electromagnetic radiation through the extragalactic medium. The medium is modelled as a large 
. number of magnetic domains, uncorrelated with one another. We obtain an analytic expression for 
the different Stokes parameters in the limit of small mixing angle. The different Stokes parameters 
are found to increase linearly with the number of domains. We also verify this result by direct 
numerical simulations. We use this formalism to estimate the effect of pseudoscalar-photon mix- 
ing on the Cosmic Microwave Background (CMB) polarization. We impose limits on the model 
^-j- . parameters by the CMB observations. We find that the currently allowed parameter range admits 
! a CMB circular polarization up to order 10~'. 
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1 Introduction 

> 

^ . Pseudoscalars with very small mass and very weak coupling with visible matter arise in many 
extensions of the standard model. The most well known example is the axion [TJ [2j [3j HJ El El [7J . 
The subject of pseudoscalar-photon mixing in background magnetic field has been studied by many 
authors [9l [lQl [TTl [T2J, [13l [13] - The mixing changes the intensity as well as the state of polarization 
of the photons. This phenomenon has also been used to search for very low mass pseudoscalars 
by laboratory experiments as well as by astrophysical observations. This search has put stringent 
limits on the mass and coupling parameters of such pseudoscalars [151 EH El HH HHJ EH 1211 [22] • 
The astrophysical and cosmological consequences have been studied extensively in the literature 
[221 [231 [2H [251 [2SI [271 EH [291 EDJ. 

The polarization of electromagnetic waves can also change due to the presence of a background 
pseudoscalar field. This phenomenon also has many astrophysical consequences [3T] 152]. 

The propagation through intergalactic medium has been studied earlier by several authors 
[2U (2S1 123 12H] • The medium is expected to be turbulent and it is reasonable to model it as a 
large number of uncorrelated magnetic domains. The magnetic field in a particular domain is 
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assumed to be uniform but points in different, random, directions as we move from one domain 
to another. Earlier studies were primarily interested in the effect of mixing on the intensity of 
the wave. In the present paper we are interested in determining the effect on all of the Stokes 
parameters. This has a wide range of astrophysical and cosmological applications in radio [25| [33] , 
CMB [22], [271 [34] and optical [25l [35l [36l [37] . Polarization is a particularly sensitive probe of the 
pseudoscalar-photon mixing [251 EE] • We assume that all the magnetic domains have equal length 
and the plasma density is assumed to be uniform throughout. A detailed analysis of the mixing, 
taking into account turbulent plasma density, is given in [TBI [25] . 

We apply our results to determine the effect on cosmic microwave background radiation. The 
mixing affects the intensity of radiation as well as the linear and circular polarization. The mixing 
also results in a distortion of the cosmic microwave background (CMB) spectrum. CMB observa- 
tions [39] can lead to stringent limits on the pseudoscalar-photon mixing, given current estimates 
of the intra-galactic magnetic fields and plasma density. 

We first set up the basic equations [38]. We consider electromagnetic waves propagating in a 
background magnetic field. The action can be written as, 



S = <Tx 



+ (1) 

Maxwell's equations excluding the effects due to gravitation are given by: 

V-E = g4> V4>-(B + B)+p; (2) 
VxE + — = 0; (3) 

VxB-^ = gjExV<f>-(B + B)^)+j; (4) 

V-B = 0. (5) 

Here Bi and Ei are the usual magnetic and electric fields and B is the background magnetic field. 
We assume that B is independent of space and time and hence set its derivatives to zero. 
The pseudoscalar field's equation of motion is 

_| _ V 2 + = —g<j)E • (B + B). (6) 

We choose the coordinate system such that the z-axis points along the direction of propagation 
and the x-axis is parallel to the transverse component of the background magnetic field B. As 
shown in Ref. [38], the longitudinal component of the background magnetic field plays a negligible 
role and can be ignored. We assume a quasi-monochromatic wave and hence the terms which are 
quadratic in the fluctuating fields give negligible contribution. We may write the field equations 
as 
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where A = E/uj and An refers to the component parallel to the transverse background magnetic 
field. The "mass matrix" or "mixing matrix" is 

M=( U l -9*Bt^\ j (8) 



where up is the plasma frequency, the pseudoscalar mass and Bt is the magnitude of the 
transverse component of B. The plasma frequency is given in terms of the electron density, n e , by 

2 4vran e 

u P = . 9 

— » 

The component of A along the y-axis does not mix with the pseudoscalar field. The matrix 
M is diagonalized by an orthogonal transformation, OMO T = Mp>, where Mp> is diagonal. We 
parameterize the orthogonal matrix, 



O 

The angle 9 is given by 



cos 9 — sin 9 
sin 9 cos 9 



tan20 = Z^£ T , (10) 
where the symbol / denotes the oscillation length, 

l = n^- (ID 

We assume that the pseudoscalar mass is very small, <C up. If the mass is much heavier than 
this, then their mixing with photons produces negligible effect for intergalactic propagation for the 
range of allowed parameters. The two eigenvalues of the matrix M, /i±, may be expressed as, 

ujl + ml . 1 



/4 = -^-^ ± - m D 2 + (2^t^) 2 . (12) 



2 Basic equations of polarization propagation 

The basic field equations, given in the Introduction, are solved by the procedure described in Ref. 
[38] . Here we simply quote the final results for propagation of the correlation functions, which are 
the density matrix elements of the polarization. The most general form of the different correlation 
functions after propagation through distance z, written out explicitly as coupled equations in 
Appendix A, can be summarized in matrix form as 

p(z) = P(z)p(0)P(z)-\ (13) 

where the matrix p(0) is defined by 

/ < A„(0)Afj(0) > <A||(0)i41(0)> <A||(O)0*(O)>\ 
p(0)= <A ± (0)A|j(0)> <A L (0)AI(0)> < A ± (0)r(0) > , (14) 
\ < 0(0) A* (0) > < 0(0)^1(0) > <0(O)0*(O)>/ 
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and the unitary matrix P(z), the solution to the field equations for a given mode u, is given in the 
case that the coordinates are chosen so that the "parallel" axis lies along the transverse component 
of the external field, Bt, by 

(1 — 7 sin 2 6* 7 cos 9 sin 9 \ 

e -«(,+A,-(^-^)'/2), q _ ^ 

7 cos 9 sin 9 1 — 7 cos 2 9 J 

In the above equations, A = A^ — A^, while 7 = (1 — e lAz ). A A and A^ are defined in terms of 
the frequency, u, and the eigenvalues, 



A A = Juj 2 - /!% - to, 



We may approximate A as 



A = Jlu 2 -fj, 2 _-u. (16) 



A = A^ - A A « - \fl + tan 2 29 w -. (17) 



The approximation, applicable to our work here, is that 3> cup, m^, and g^Bx- We also find that 
in our case tan 29 C 1, as we elaborate below. For later reference, we give the leading term in the 
expansion of the phase in the 22 element of Eq. (fT5l) . which we will denote 8z, in terms of sin#, of 
order 10 -5 in our application: 

iSz = i{u + A A - {u 2 - co 2 P ) 1/2 )z w -i~ sin 2 9, (18) 

where up u for n e ~ 10~ 8 cm" 3 , as commonly estimated [23], and u ~ 50 GHz for the CMB. 
The consequent expansion of the 22 element of P{z) reads 

e~ i5z w 1 + iy sin 2 ^. (19) 

In the model that we consider, the medium has a uniform value of up, namely uniform electron 
density, and the cluster diameters are all the same, with value z, over which the direction and 
strength of the magnetic field is assumed to be uniform. Each cluster has a random orientation 
with respect to some "external" coordinate system, where the propagation direction is z. Eq. (fl3l) 
may be taken as the propagation from the initially prescribed polarization correlation matrix p(0) 
through a distance z, in a uniform field of strength Bp that is oriented along the "parallel" axis. 
We take the magnetic field in the first cluster to be aligned at an angle f3\ to the fixed external 
coordinate system. Therefore we rotate the initial state by an angle f3\ with rotation matrix R(Pi) 
to align the "parallel" axis along the magnetic field direction in the first cluster. After a distance z, 
a new zone (cluster) is entered, where the field has the same strength but a random orientation with 
respect to the first. At this point, after the propagation through the first cluster, the polarization 
state has changed. The new state is still expressed in terms of the rotated coordinate system, of 
course. So we rotate back by angle (3\ to return to the frame of the external coordinate system. 
Now the state is rotated by an angle fa with rotation matrix R(fa) to align the new "parallel" 
axis along the magnetic field direction in the second cluster. The propagation through distance 
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z in the second cluster results in a new polarization density matrix, expressed now in terms of 
the new rotated coordinates. Rotate back to the original coordinate system, then forward to the 
new (third) cluster's magnetic field direction. This process continues until the propagation wave 
reaches "Earth" , after, say, n clusters have been traversed. The propagation then amounts to the 
following matrix product: 

Pn (z) = R-\(3 n )P(z)R((3 n )R-\^ 1 )P(z)R((3 n ^)R~\(3 n . 2 )..M- 1 ((3 1 )P(z)R((3 1 ) 

x P ^)R-\p l )p-\z)R(p l )R-\p 2 )p-\z)R(p 2 )...R-\(3 n )p-\z)R(f3 n ) 1 (20) 

where the rotation matrix R{j3 m ) acts only on the two dimensional space transverse to the propa- 
gation direction, and it reads 



R(Pm) = -sin/3 m cosP m . (21) 




Further simplifications occur, because i?(/3 m )i? _1 (/? m _i) = R(/3 m — /3 m -i), since R~ 1 {(3)=R{—[3). 
The reference axis for (3 is arbitrary, since only the relative angles between the field directions 
in adjoining clusters are relevant. The difference between two random angles is obviously a 
random angle, so the structure of Eq. (1201) amounts to the unitary transformation p n (z) = 
U(z,n)p(0)U~ l (z,n), where U(z,n) is the product of random rotations of P(z). Next, since the 
parameter 9 is very small (~ 10~ 5 in our case), we develop Eq. (120]) to the leading order in 9. 

To expand the propagation matrix in powers of sin. 9, write it as the sum of a (nearly) unit 
matrix and a matrix proportional to sin 9: 

P(z) =l(z) + sin 9 Viz, 9). (22) 

Suppressing the irrelevant over all phase in Eq. ( fl5l) . the matrices 1{z) and V(z, 9) are given by 




X(z) =(01 I (23) 
and 

(—7 sin #07 cos 9 \ 
zfsintf , (24) 
7 cos 9 7 sin 9 J 

and 7 is defined as in Eq. (fT5|) . Using the fact that T(z) and R(0) commute and X(z)^ 1 and 
R{(3)~ 1 = X(—z) and R(—(3), respectively, one finds that to leading order in sin^ V(z,9), 

pl\Z) = l(zTp(0)l(-zT 

n 

+ sin 9 ^ [2{z) n ' j R{-f3 3 )V{z, 9)R{(3 j )l(zy- 1 p(0)l(-z) n + h.c] , (25) 
i=i 

where Z = nz is the total distance travelled. The form of p n ~ pH\ Eq. (|25p . a leading zeroth 
order term and a sum of pairs of first order terms, each corresponding to non trivial propagation 
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in only one of the n cells, is what one expects to leading order. The above treatment of the model 
is perfectly general. 

We next assume that the electromagnetic wave is initially unpolarized. Furthermore we assume 
that the initial correlators, 

< 0(O)AJj(O) > = < <f>(0)A* ± (0) > = < 0(O)0*(O) > = . (26) 

Hence the initial density matrix, 

p(0)=diag(l,l,0). (27) 

In this case the order sin# terms, as defined in Eqs. (1221) and ft2M . are actually of order sin 2 9 in 
the terms p(z)n, p(z)i2, p(^)2i and p{z)ti that affect the polarization parameters at first order in 
sin 9 V(z,9). The j th term of the sum in the expression for p^ reads, for the special case where 
p(0) = (1,1,0), 

/-2sin 2 ^cos 2 /? i (l - cosAz) - sin 2 9 sin 2(3,(1 - cosAz) 7 * e - iAz ( n -J') cos sin(2#)/2 
p W( z ) i= - sin 2 9 sin 2/3,(1 - cos Az) -2 sin 2 9 sin 2 ^(1 - cos Az) 7 * e -« A *("-i) sin/3,- sin(2#)/2 
\ 7 e iAz (""J) cos Pjsm(20)/2 7e iA ^-i) s i n ^ sin(2fl)/2 

(28) 

At second order in sin^ V(z, 9), the mixing between A|| and 0, and Aj_ and induces order sin 2 9 
terms, so these must be included as well, and we deal with them next. 

There are 2n(2n — l)/2 ways to choose two propagation factors V(z, 9) from the 2n positions in 
the general expression Eq. (1201) . The general second order expression is an extension of the form 
shown in Eq. (1251) . Denoting it p n , we find 

n k— 1 

P V\Z) = sin 2 9[ E P(^) n ~^(-A)^(^ 9)l(z) k - j ~ 1 R(f3 k - fr)V(z, 9)R{^)l{ Z y- 1 

k=2 j=l 

x p(0)X(-z) n + h.c] 

n n 

+ E E^r^i-^)^ 9)R(p ] )i( Z y- i p(o) 

k=l j=l 

x I(-z) k - 1 R(~(3 k )V(-z,9)R(p k )I(-z) n - k ]y (29) 

Using Eq. (|23|) . X(z) m is simply diag(l,l,e iAm2; ), so i?(/3) and X(z) commute. Thus their orders 
in Eqs. (1251) and (|29|) are immaterial and can be chosen for convenience. Only the propagation 
matrix V(z, 9) and, in principle, the initial conditions matrix, p(0), do not generally commute with 
each other and with the R(0) and 1(z) matrix. In general their positions in the matrix product 
are fixed, though for special cases such as ours, which has the unpolarized initial condition with 
zero value for the field initially, p(0) = diag(l,l,0), the matrices can be collapsed and simplified 
even further, as illustrated in Eq. (1281) and below in Eq. (130|) . 

When p(0) = diag(l,l,0), the p(z) elements relevant to our study of spontaneous polarization at 
leading order are of the form 7 2 sin 2 9 x products of cosines and sines of rotation angles. Corrections 
to the leading terms are themselves of order sin 2 9 w 10~ 10 for the parameters we choose for our 
model. Assembling all of the contributions of terms of second order in sin^ V(z,9), we find that 
the j-, k- th term of the sum in Eq. (|29|) . yields the following elements of the 2x2 submatrix 
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relevant to polarization: 



Pn (z,j, k)n = —4 sin 2 6(1 — cos Az) cos[A^(A; — j)] cos fa cos fa 

p [2] (z,j,k) 12 = -2sin 2 6>(l - cos Az) cos[Az(k - j)](sin(fa + fa) -isin(fa - fa) t&n[Az(k - j))) 

p [ n ] (z,j, k) 21 = -2 sin 2 0(1 - cos Az) cos[Az(k - j)](sin(fa + fa) + % sin(fa - fa) tan[Az(A; - j)}) 

p$(z,j,k) 22 = — 4 sin 2 6(1 — cos Az) cos[Az(k — j)] sin fa sin fa. (30) 

It turns out that all of these contributions come from the terms in the first double sum in Eq. 
since the zero in the 33 entry of p(0) prevents the order one elements 13, 23, 31 and 32 in 



j from communicating with those in p L n(z)k to form order one elements in the 11, 22, 12 and 



21 product in the second double sum term in Eq. (1291) . 

At this point, combining each term with its hermitian conjugate, we have displayed n(n — 
l)/2 + n = n 2 /2 + n/2 terms at order sin 2 6. It is this second order term in sin# V(z,6) that 
dominates the CMB photon's random walk through the clusters. Inspection of Eq. ( |30l) reveals 
that the polarization parameters Q, U and V are, on average, all the same order of magnitude. 
Our analysis leads to the conclusion that the number of contributing "steps" in the random walk 
represented by the propagation process is proportional to n 2 , so the effective displacement will 
be proportional to n, with the spontaneous polarization at a step being of order 10~ 10 . The 
dependence of the intensity on the number of domains has been discussed earlier [241 [26l [27J [28] . 
However the linear dependence of polarization observables Q, U and V on n is a new result. For the 
transport through a thousand cells, or clusters, we expect a typical net displacement (polarization) 
of order 10~ 7 , given our parameter assumptions. 

We can compute reduced Stokes parameters using the correlation functions of A||(2) and (f>(z), 

I = < A\\(z)A\\(z) > + < A ± (z)A* ± (z) >= p u + p 22 , (31) 

Q = <A l] (z)Al(z)>-<A x (z)A* ± (z)>=p 11 -p 22 , (32) 

U = <A ± (z)Al(z)> + <A ll (z)A* L (z)>= Pl2 + p 21 , (33) 

V = i(-<A ± (z)AUz)> + <A ll (z)A* ± (z)>) = i(-p 12 +p 21 ). (34) 



The linear polarization angle is given by, 



tan2^ = ^, (35) 



and the degree of polarization is given by, 

y/Q 2 + W + V 2 

P=- j • (36) 

With this background analysis of the model, we turn to applications to CMB temperature and 
polarization. 



3 CMB Temperature and Polarization 

In this section we consider the effect of pseudoscalar mixing on the propagation of CMB photons. 
The mixing affects both the CMB temperature and polarization. We can impose a limit on the 
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pseudoscalar-photon coupling by requiring that the CMB is distorted to less than one part in 10 5 , 
which is the amplitude of CMB temperature fluctuations. Furthermore we also demand that the 
mixing does not generate a CMB polarization larger than 10 -6 . Finally since the WMAP has not 
observed B modes at the level of 10~ 6 , we require that these are at least an order of magnitude 
smaller. The limit can be imposed only if we assume a value for the background magnetic field 
and plasma density. The magnetic field, in particular, is very poorly known. Furthermore the 
results can change significantly if the magnetic field and plasma density changes along the path 
as well as if there is a background flux of pseudoscalars. We set the incident pseudoscalar flux to 
zero. Furthermore we assume a uniform background field and plasma density within a particular 
domain. The result is also found to depend significantly on the domain size that we assume and 
hence on the details of how we model in the intergalactic medium. Due to the uncertainty in the 
magnitude of the magnetic field it is best to impose a limit on the product g^Bx- This is possible 
since the magnetic field always occurs in the form of this product. The model dependence implies 
that the limit also depends on some parameters of the model, such as the domain size and/or the 
total number of domains. We shall address this issue later in this section. 

The most stringent limits on the coupling are obtained from astrophysics [15]. The current 
astrophysical limit is < 6 x 1CT 11 GeV -1 . The typical upper limit value of the intergalactic 
magnetic field is about 10~ 9 G [ID] and the electron number density n e ~ 1CT 8 cm -3 [21]. We 
point out that the supercluster magnetic field may be significantly larger, of order 10~ 7 G [U]. In 
broad outline, we model the intergalactic medium with several magnetic domains, each of size z 
of order a few Mpc. The total propagation distance is of the order of a few Gpc. We shall assume 
that the magnetic field in different domains is uncorrelated. 

Choosing specific values, we calculate the effect on CMB temperature and polarization observ- 
ables as CMB radiation propagates through n = 1000 magnetic domains, each of size z — 1 Mpc. 
The relevant physical quantities are: 

1. Magnetic field, B T = 10' 9 G, 

2. Electron number density, n e = 10~ 8 cm -3 , 

3. Frequency of radiation (CMB), v = 50 GHz, hence the oscillation length parameter, I = 200 
pc. 

We set the coupling g$ — 6 x 10~ n GeV -1 . We numerically compute the Stokes parameters after 
propagating through n clusters by using the equations given in Appendix A. After propagating 
through each cluster we rotate the electromagnetic wave vector in order to account for the change 
in the direction of the transverse magnetic field from one cluster to the next. For this purpose we 
need to transform the correlators in a fixed reference frame, at the begining of each cluster, to the 
reference frame with x-axis parallel to the transverse magnetic field in the cluster. The relevant 
formulas for the transformation of the correlators are given explicitly in Appendix B. We point 
out that the initial density matrix is given by Eq. f[2"T|) . 

In Fig. [I] we show histograms of 1000 simulations of Iq — I, Q, U and V Stokes parameters 
after propagating over 1000 clusters. All the parameters are normalized by the initial intensity Jo- 
Here we have assumed that initially the wave is unpolarized. If initially the wave has a non-zero 
value of the Stokes Q parameter, our results remain essentially unchanged with the mean value of 
the Q parameter shifted by its initial value. The effect on the U and V polarization is negligible. 
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In CMB studies one deals with the coordinate independent E and B modes. In the present 
case we work directly with the coordinate dependent Q and U Stokes parameters. This is fine as 
long as we are working in a particular direction in the sky. We align the coordinate system such 
that in the absence of pseudoscalar-photon mixing, only the Q polarization is non-zero. Then the 
U polarization generated by the mixing in the chosen frame provides an estimate of the B modes. 
Hence the polarization generated by the pseudoscalar-photon mixing must be less than the current 
limit on the B modes. 

In Fig. [1] the intensity is found to deviate from its inital values by less than 1 part in 10 6 
due to mixing for the chosen parameters. Hence the effect of mixing on temperature anisotropy 
is negligible. Also the Stokes parameter Q remains very close to its initial value of zero. We see 
that the predictions are not in conflict with the current observations. We predict relatively large 
fluctuations in the U parameter and the circular polarization V. The sample mean values of U and 
V are found to be approximately zero, typically of the order of 10~ 9 . The corresponding standard 
deviations are found to be 2.0 x for both of the parameters. Hence we find that the choice of 
parameters already predict the fluctuations in the U polarization close to its limiting value. We, 
therefore, impose a limit on the product g^Br < 0.2 Mpc -1 due to CMB polarization. 

In Fig. [2] we show the fluctuations in Q, U and V as a function of distance for a particular 
realization of the model. The plot for this individual case shows similar magnitudes induced for 
all the three Stokes parameters. This is true for the entire distribution of realizations as well, as 
shown in Fig. [TJ 

With results of the full simulation in hand, we can check how the results depend on the number 
of clusters. In Fig. Owe plot the dependence of mean (absolute value) and standard deviation of 
U on the number of clusters. The cluster size is kept fixed in this calculation. As we anticipate 
from our discussion in Sec. 2, we find that both the mean and standard deviation increase linearly 
with the number of clusters. The standard deviation shows very little digression from the linear 
plot, whereas the mean shows relatively large fluctuations. We point out that this linear behavior 
is expected as long as the Stokes parameters are much smaller compared to unity. If the number 
of clusters or the magnetic field in each cluster becomes too large then the linear dependence will 
break down. 

We next determine how our predictions depend on the assumed model of intergalactic space. In 
Fig. H]we show how the mean and standard deviation of the U parameter change as we reduce the 
size of each cluster and proportionately increase the total number of clusters such that the total 
distance of propagation remains unchanged. We find that the U Stokes parameter increases as the 
cluster size is reduced. This is easily understood. The oscillation length I is much smaller than the 
size of each cluster. Hence, statistically, the contribution per cluster remains almost unchanged 
as we decrease the size z of the cluster as long as z 3> I. This is clear from the analytical result, 
Eq. fl5U|) . obtained in Sec. 2. The result depends on z only through the functions such as cos Az, 
ta,n[Az(k — j)]. Since the arguments of the trigonometric functions are much larger than unity, 
adding a large number of such terms yields a negligible or fluctuating dependence on z. However 
the total number of clusters, n, keeps increasing, leading to a linear increase with n in the total 
contribution, as explained in Sec. 2. Due to this model dependence we may express the limit more 
reliably as 



We expect this to be valid as long as the cluster size, z, is much larger than 1/A. Although we 




(37) 
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Figure 1: The distribution, f(x) (xlCT 6 ), of the normalized Stokes J, Q, U and V parameters for 
the CMB. For / we only show the difference Jo — I, where I is the initial intensity. All of the 
parameters have been normalized by dividing by I . Each histogram contains the results of 1000 
simulations. 
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Figure 2: The distance dependence of the normalized Q, U and V Stokes parameters for a particular 
random realization. 
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Figure 3: The mean and standard deviation of U/Iq as a function of the number of magnetic 
clusters n. For each n we use a total of 1000 simulations. Both the mean (lower curve) and 
standard deviation (upper curve) show linear dependence on n. 
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Figure 4: The mean and standard deviation of U/Iq as a function of the number of magnetic 
clusters n. Here the total distance of propagation is kept fixed. For each n we use a total of 1000 
simulations. Both the mean (lower curve) and standard deviation (upper curve) show roughly 
linear dependence on n. 

have obtained this assuming uniform cluster size, and magnetic field strength equal in all clusters, 
the result is approximately valid even if the magnetic field and cluster size fluctuates from cluster 
to cluster. In this case we interpret the magnetic field strength in Eq. (1371) as the mean magnetic 
field over all the clusters. We have explicitly verified this numerically by allowing fluctuations 
in the magnetic field and cluster size with standard deviation of the order of these parameters. 
However the result may deviate significantly if these parameters show much larger fluctuations. 

What part of the density matrix is driving the linear growth with the number of clusters n? 
From our discussion in Sec. 2, we see that the linear growth comes from the second order terms 
which involve the correlators of </> with An and Aj_. We can verify this numerically by setting all 
such correlators to zero at each step. In Fig. [5] we show the result of this numerical experiment. 
The consequence is a growth in the Stokes parameter U proportional to y/n. This confirms our 
expectation in Sec. 2 that the linear growth is driven by the correlators < (p(z)A^(z) > and 
< ®{z)A* (z) >. 

The values of the U and V Stokes parameters are found to be relatively large. Such large values 
of the U and V parameters are somewhat surprising and it is useful to get an independent estimate 
to verify the numerical results. For the parameters chosen the mixing angle \6\ ~ q^BtI/Z ~ 10~ 5 . 
Furthermore the oscillation length I « 10~ 4 Mpc <C z. Hence we expect that in one domain the U 
polarization would be of order 8 2 ~ 10~ 10 . As we have seen above the effect is linearly proportional 
to the number of clusters. Hence for 1000 clusters we multiply this number by 1000 to get 10~ 7 , 
in agreement with the value found by direct numerical computation. 
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Figure 5: The mean and standard deviation of U/Iq as a function of the number of magnetic 
clusters n. Here all the correlations of with fields other than have been artificially set to zero 
at the beginning of each cluster. For each n we use a total of 1000 simulations. Both the mean 
(lower curve) and standard deviation (upper curve) show linear dependence on y/n. 

4 Summary and discussion 

Proposed as an extension to the standard model in order to solve the strong-CP problem in 
quantum chromodynamics (QCD), axions today still remain undetected. Direct and indirect (such 
as ours) methods of detection are being used to search for axion-like particles. Future missions 
such as Planck will improve the sensitivity to the polarization of the CMB, which is a crucial next 
step in pinning down the origin of anisotropies in the CMB. A detection of circular polarization 
would certainly strengthen the case for existence of ultra-light pseudoscalars (such as axions). 

In this paper we have analyzed the mixing of pseudoscalars and photons in order to explain 
the polarization of the cosmic microwave background radiation. We have treated the intergalactic 
region as a collection of magnetic domains (each of size 1 Mpc), and have shown, both theoretically 
and numerically, that the polarization of radiation increases linearly with the number of clusters 
(domains). Within the experimentally allowed range of g^Bx values, our model predicts a B mode 
polarization and circular polarization up to order 10~ 7 . 

The model that we have used assumes that magnetic fields in adjoining clusters are completely 
uncorrelated. Although this is the most reasonable assumption, it would be interesting to see how 
things change if the fields are instead correlated. We might, for example, be able to explain the 
reported observation of large scale alignment of the optical polarization of quasars [53] . We intend 
to address this issue in a future publication. 
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Appendix A 

Spelled out in components, the propagation equation of the density matrix p(z) amounts to six 
independent equations: 



< A\\(z)Al(z) > = - < A||(0)Afj(0) > 1 + cos 2 26 + sin 2 26 cos[,2(A^ — A^)] 



+ 
+ 



< A ± (z)A* ± (z) > 
<A ± (z)Al(z)> 

< <j>(z)4>*(z) > 



< 0(O)0*(O) > sin 2 26 - sin 2 26 cos[z(A - A A )} 
]- < 0(0)^(0) > [sin 26 cos 26 - sin 26 cos 26 cos[^(A^ - A A )} 

i sin 26 sin[^(A - A A )]J + c.c.j 
< A ± (0)A* ± (0) > 



= < A ± (0)Afj(0) > 
+ < A ± (O)0*(O) > 



cos" 6 e 1 Fz + sin" 6 e 1 



Gz 



sin 6 cos 6 (e 



i „i Gz 



(38) 
(39) 

(40) 



2 
1 



< A||(0)Afj(0) > sin 2 26 - sin 2 26 cos[z(A - A A ) 



+ ^ < 0(O)0*(O) > 1 + cos 2 20 + sin 2 20 cos[^(A^ - A A )} 



+ 



< (f>(0) ATM > -sin 26 cos 26 + sin 2# cos 26 cos[^(A - A A ) 



+ i sin 26 sin[^(A - A A )\ 



c.c. 



(41) 



< (j>(z)A\(z) > = - < A||(0)Afj(0) > sin 2$ cos 2$ — sin 26 cos 2$ cos^A^ — A^)] 
— i sin 2$ sin[z(A0 — A^)] 

+ - < 0(O)0*(O) > —sin 26 cos 26 + sin 2$ cos 26 cos^A^ — A^)] 
2 

+ i sin 26 sin[^(A - A A )} 

+ - < 0(0)^(0) > [sin 2 26 + (l + cos 2 26) cos[z(A - A A )} 
2 

+ 2i cos 26 sin[z(A^ - A A )} 

+ - < A||(O)0*(O) > [sin 2 26 - sin 2 26 cos[^(A - A A )]1 
= <A||(0)A1(0)> ^ 

+ < 0(0)^1(0) > [sin 2 6 e" 1 Fz + cos 2 6 o" 



sin 6 cos (e" 1 Fz - e~ 1 Gz ) 

i Fz , „ 2 /i - i Gz 



(42) 



(43) 
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where 

F » ^l-"D (44) 
G » ^--^p) (45) 
and the remaining symbols are defined in Sec. 1 and 2 in the text. 

Appendix B 

Let Ai(n, 0) and A 2 (n, 0) be the two components of the electromagnetic field in a fixed reference 
frame at the begining of the n th cluster. For each cluster we define a local reference frame with 
x-axis aligned parallel to the transverse component of the background magnetic field. In the 
local frame corresponding to the n th cluster we denote the two electromagnetic field components 
as 0) and A±(n, 0). We assume that initially, i.e. at n = 0, z = 0, only the correlators 

< A 1 (0,0)AI(0,0) > and < A 2 (0, 0)^(0, 0) > are nonzero. After propagation through the first 
cluster all the six correlation functions are likely to be non-zero. After transforming into the local 
frame, the initial correlators for the n th cluster are, 

< A\\(n,0)A*\(n,0) > = cos 2 a < A 1 (n,0)A* 1 (n,0) > + sin 2 a < A 2 (n,0)A 2 (n,0) > 

+ sin a cos a (< A 2 {n,0)A* 1 {0) > + < A 1 (n,0)A* 2 (0) >), (46) 

< A±(n, 0)A* ± (n, 0) > = sin 2 a < A^n, 0)Al(n, 0) > + cos 2 a < A 2 (n, 0)A* 2 (n, 0) > 

- sin a cos a (< A 2 (n,0)A* 1 (n,0) > + < A 1 (n,0)A 2 (n,0) >), (47) 

< A ± (n, 0)A* l (n, 0) > = cos 2 a < A 2 (n, 0)A*(n, 0) > - sin 2 a < A^n, 0)A*(n, 0) > 

+ sin a cos a (< A 2 (n,0)A* 2 (n,0) > - < A 1 (n,0)A* 1 (n,0) >), (48) 

< 0(n,O)Aj|(n,O) > = cos a < (f)(n,0)A* 1 (n,0) > + sin a < (f)(n,0)A 2 (n,0) >, (49) 

< (f)(n, 0)A* ± (n, 0) > = - sin a < <j>{n, 0)A{{n, 0) > + cos a < <f>(n, 0)A*(n, 0) >, (50) 

and < 0(n, O)0*(n, 0) >. 
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